Environmentally induced corrections to the geometric phase in a two-level system 
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We calculate the geometric phase for different open systems (spin-boson and spin-spin models). 
We study not only how they are corrected by the presence of the different type of environments 
but also discuss the appearence of decoherence effects. These should be taken into account when 
planning experimental setups to study the geometric phase in the nonunitary regime. We propose 
a model with slow decoherence rate in which the geometric phase is still modified and might be 
measured. 



Since the work of Berry [l|, the notion of geometric 
phases has been shown to have important consequences 
for quantum systems. Berry demonstrated that closed 
quantum systems could acquire phases that are geometric 
in nature. He showed that, besides the usual dynamical 
phase, an additional phase related to the geometry of the 
space state is generated during an adiabatic evolution. 

The existence of such a phase is also true for open 
quantum systems. In particular, when a static potential 
is exerted on the main system, the wave function of this 
system acquires a phase and hence the interference term 
appears multiplied by a phase factor e lv . In an interfer- 
ence experiment, its effect on the pattern of the system is 
related to the phase's statistical character, particularly, 
in situations where the potential is not static. Yet more 
importantly, any source of stochastic noise would create 
a decaying coefficient, usually called decoherence factor 
F. For a general case, the phase if is described by means 
of a distribution function 0, H[ ■ No matter how weak the 
coupling that prevents the system from being isolated, 
the evolution of an open quantum system is plagued by 
nonunitary features like decoherence and dissipation. De- 
coherence, in particular, is a quantum effect whereby the 
system loses its ability to exhibit coherent behaviour and 
appears as soon as the partial waves of the main system 
shift the environment into states orthogonal to each other 
Nowadays, decoherence stands as a serious obstacle 
in quantum information processing. 

The geometric phase (GP) for a mixed state under 
nonunitary evolution has been defined by Tong eia/.[f| 
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where e k {t) are the eigenvalues and \^k) the eigenstates 
of the reduced density matrix p r (obtained after tracing 
over the reservoir degrees of freedom). In the last defini- 
tion, t denotes a time after the total system completes 
a cyclic evolution when it is isolated from the environ- 
ment. Taking the effect of the environment into account, 
the system no longer undergoes a cyclic evolution. How- 
ever, we will consider a quasicyclic path V : t e [0, r] 
with r = 27r/0 (f2 the system's frequency). When the 



system is open, the original GP, i.e. the one that would 
have been obtained if the system had been closed 
is modified. That means, in a general case, the phase is 
$ = <& u + (5$, where <5$ depends on the kind of environ- 
ment coupled to the main system @. 

It is expected that GPs can be only observed in in- 
terference experiments carried out in a time scale slow 
enough to ignore nonadiabatic corrections, but rapid 
enough to avoid the destruction of the interference pat- 
tern by decoherence Q- So far, there has been no ex- 
perimental observation of GPs for mixed states under 
nonunitary evolutions. The purpose of this short article 
is to study how GPs are affected by decoherence in dif- 
ferent physical scenarios. The decoherence time results 
very important when trying to measure the GPs since 
for times longer than the former the GPs, literally, dis- 
appear. In this framework, we shall compute the GP for 
different models using the kincmatical approach to the 
GP given by Eq.fl]), and compare the results therein ob- 
tained. We shall start by reviewing some of our previuos 
results [H, and then we shall present further results con- 
cerning the environmentally induced corrections to the 
GP (S$) in realistic (even experimentally feasible) mod- 
els. 

Purely Decohering Solvable Spin-Boson Model. In this 
section, we shall review the basic results for an open 
quantum system by presenting a model which is simple 
enough to be solved analytically Q. In spite of its sim- 
plicity, this model captures many of the elements of deco- 
herence theories and sheds some insight into the modifi- 
cation of the GPs due to the presence of the environment. 
This model has been used by many authors to model de- 
coherence in quantum computers[9( and, in particular, it 
is extremely relevant to the proposal for observing GPs 
in a superconducting nanocircuit [Iol |. The Hamiltonian 
that describes the complete evolution of the two-state 
system interacting with the external environment is: 
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where the environment is described as a set of harmonic 
oscillators with a linear coupling in the oscillator coordi- 
nate. The interaction between the two-state system and 
the environment is entirely represented by a Hamiltonian 
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in which the coupling is only through a z . In this partic- 
ular case, [<r z ,Hi nt ] = and the corresponding master 
equation is much simplified, with no frequency renormal- 
ization and dissipation effects. In other words, the model 
describes a purely decohering mechanism, solely contain- 
ing the diffusion term T>(t) whose master equation, after 
tracing out the environmental degrees of freedom, is given 
by (with H = 1) 
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where V(s) = J Q S ds' J °° dojl(uj) coth ^ 2fc ^ T ^ 
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s')), and I(u>) is the spectral density of the environment, 
usually, I(lu) ~ w" up to some frequency A that may be 
large compared to f2. In particular, the case with n = 1 
is the "ohmic" environment. 

Then, it is easy to check that p Tol (t) = 
e - lf2 *-^(t)p roi (0) is the solution for the off-diagonal 
terms (while the populations remain constant), where 
A(t) = LdsT>(t). In the following, we shall call 
F = cxp(— A(t)) the decoherence factor. 

Hence, the GP for an initial pure state of the form 
*(0)) = cos#o/2|e) +sm8 /2\g), related to a quasicyclic 
path V : t e [0, r] up to first order in the dissipative 
constant (70 oc A 2 ,) is[8j 

$sb ~ 7r(l - cos# ) 
+ O(7o 2 )- 
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In the right side of last expression, we have performed 
a serial expansion in terms of 70. The first term corre- 
sponds to the unitary phase Consequently, we see 
that the unitary GP is corrected by a term which de- 
pends directly on the kind of environment present Jj|. 
For example, for an ohmic environment in the limit of 
high temperature <5<1>c?b = 71-2 {Io/^I^bT sin 2 9q, while 
the same environment at zero temperature modifies the 
unitary phase as £$sb = f 7o( — l + log(27rA/JT)) sin 2 9 . 
These results can be compared with those in [j, [ll[ . In 
those cases, the correction due to the environment is also 
proportional to (70/^) sin 2 (#o) (mainly due to the simpli- 
fied decoherence factor F = exp(— jot)). However, in our 
model, these corrections enclose the main characteristic 
of the model of bath we are taking into account, which 
allows to evaluate the decoherence time scale properly. 

In the case of having a bosonic environment, composed 
by an infinite set of harmonic oscillators, it is not diffi- 
cult to evaluate the decoherence time scale. This scale 
should be compared with the time r = 27t/tt at which 
one expect to measure the GP. In the case of an ohmic 
bath in the high temperature limit, the decoherece time 
is to = l/(7o7rfcsT), which is really a very short time 
scale compared with r. In the zero temperature case, the 
decoherence time scales as ~ e 1 / 7o /A which, indeed, 
can be very large in the case of underdamped environ- 
ments. In conclusion, one could expect that the GP can 



be only detected at very low temperature when the atom 
is mainly coupled to a bosonic fieldQ| . 

Spin-Spin Model. We shall study another simple solv- 
able model in which the size of the environment has a 
relevant role. Consider a two-level system coupled to n 
other two-level systems fl2j. Our main subsystem (one 
qubit) interacts with the rest of the environmental spins 
by a bilinear interaction described by the interaction 
hamiltonian 



7T N 
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(5) 



fc=2 



where the system qubit is denoted by the superscript 
"1" . This coupling is also a purely phase damping mech- 
anism, as in the spin-boson model mentioned above. 
Given a factorizable initial state of the form |$(0))i = 
[a|0)i + 6|l)i] Hk =2 ( a k\0)k+(3k\^)k), the interacction en- 
tangles the state of the system with the environment. 
This means that after the interaction, both system and 
environment states are not longer factorizable. Similarly 
to the spin-boson model, the density matrix will have 
constant populations (since [a z , Htot] = 0) and the off- 
diagonal terms will be multiplied by a decoherence factor, 
as Pq X = ab*z(t) where 

JV 

z{t) = \\ [cos(7r Jifct) + sin(TrJifei)], (6) 

fe=2 

where <p± = |aj,| ± \(3k\- Note that z(t) depends on the 
initial conditions of the environment only through the 
probabilities of finding the system in the eigenstates of 
the interaction Hamiltonian \au\, \(3k\ [13] • I n this case, 
z(t) plays the role of the decoherence factor F since con- 
tains the information related to the tracing out of the 
spin environment degrees of freedom. In particular, the 
magnitude of z(t) determines the damping of the phase 
information originally contained in poi(O). In particu- 
lar, when \z(t)\ — > 0, the nonunitary evolution and the 
irreversibility of the process are evident. However, infor- 
mation can be in principle recoverable for a finite system 
since \z(t)\ is at worst quasiperiodic[T3]- The effective- 
ness of the decoherence mechanism is determined by the 
dimension of the environment. However, in any case, if 
z(t) is a complex function, it implies a phase shift and an 
attenuation of the interference fringes, i.e. a dephasing or 
decoherent process. In principle, the correction induced 
on the GP is the same as in Eq.(f4]), just replacing F(t) 
by z(t). 

Let's take for example the particular case when the en- 
vironment is composed of only one spin (k = 2 in Eq.([5])). 
For the same initial state mentioned above, and consid- 
ering |afc| = |/3fc|, we obtain z(t) = cos(7rJi) (where we 
set J = J12). In this case, z(t) is real and then, its only 
contribution is to the phase shift of the system, while 
one spin environment is not effective inducing decoher- 
ence on the system. Nevertheless, we will show that this 
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factor induces a correction to the GP which is quadratic 
in the coupling strength with the environment. In such 
a case, if one performs a serial expansion in powers of 
the coupling constant J, one obtains that the modifica- 
tion to the unitary phase is at second order. Thus, the 
correction to the unitary GP is given by 

4tt 4 

M> z ss «— J 2 sin 2 o . (7) 

This simple result shows that correction to the unitary 
GP induced by the presence of this environment can be, 
in principle, detected in an interference experiment, with- 
out the constraint imposed by the dccohcrcncc time scale. 
At zero-order, the unitary GP is the same as in Eq.Q 
<S> U =tt(1 -cos0 o )- 

Hierarchical Qubit-Qubit Decoherence Model. Herein, 
we shall compute the GP's correction for a model very 
similar to the above described spin-spin one. This sce- 
nario has the particular feature that it can be imple- 
mented to simulate quantum decoherence fl3l ]. In this 
case, the environment is also limited to only one spin 
(qubit). However, through the strategy of randomly re- 
dressing the phase of the environment qubits during the 
interaction with the system, it is possible to simulate a 
much larger environment. Therefore, the result must be 
averaged over many realizations of this evolution. The 
dimension of the Hilbcrt space can not be larger than 
N 2 , where N is the dimension of the local main sys- 
tem. To remove the information from the finite quantum 
environment, a classical stochastic field is included. Ba- 
sically, the technique consists of applying classical kicks 
to the environment qubits, and then averaging over the 
realizations of this stochastic noise. This has the effect 
of scrambling the system information after it has been 
stored in the quantum environment through the coupling 
interaction. 

We shall consider the evolution of this system subject 
to a sequence of kicks that only affect the environment 
qubit. Every kick is generated by a transverse mag- 
netic field whose effect is to rotate the environment qubit 
around the y axis by an angle e included randomly in the 
interval (—a, a). In this case, the reduced density matrix 
is similar to the above models, but for a different decoher- 
ence factor F. The off-diagonal terms are p n . = ab*fij, 
where carries all the information about the effect of 
the environment qubit on the system qubit. It is obtained 
after tracing out the environment degree of freedom and 
averaging over the many realizations of the external mag- 
netic field [H|]. In the case that there arc no kicks, i.e. 
a = 0, and /12 = cos(-7r Ji) — ip z sin(7r Ji), which agrees 
with the spin-spin model described above (p z is the ini- 
tial polarization of the environment qubit). In this case 
there is no decoherence and the GP-correction is given 
by Eq.0. The decoherence factor is independent of the 



kicking rate (no kicks in this limit), and the system qubit 
rotates independently of the environment qubit. 

If one allows a complete ramdomization, i.e. the kick 
angles tj may vary over the entire interval between 
and 27T, the decoherence factor can be approximated, in 

the limit of faster kicks, by[l3| /oi(r, t) tts e~ ™ 2'r ' — 

7 7T 2 J 2 t 

ip z sin(5^-)e sr - , where T is the kick rate. Using this 
expression, one can evaluate the correction induced on 
the GP (<5<5gs) (for the particular case p z = 0) as 

^Is«2f^ 2 sin 2 0o. (8) 

In this situation, the decoherence time is given by to = 
2T/(J 2 ) which is larger than r, making the decoherence 
process negligible if trying to measure these corrections 
to the GP. 

Finally, we shall consider the case of small angles, since 
it is the regime used by simulations and also for deco- 
herence experiments (usually with a = it/ 20). In such 
a case, it is possible to estimate the decoherence factor 
a s /12 = e~ rt£ (l + e/2)[cos(7rJf) — ip z sin(7rJi)], where 
e = 2/3a 2 is a small number (e w 0.016 for the given 
experimentally accesible value of a mentioned above). 
This decoherence factor determines a very large dephas- 
ing scale: to = l/(re). In this case, we can also evaluate 
the environmentally induced correction to the GP (up to 
sencond order in the coupling with the environment and 
also for small e, and p z = 0) <5$|g as 





2 * 4 j2 
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This correction to the GP has a term independent of 
the coupling constant with the environment J, which in 
this limit is linear with e, the small angle that is rotated 
due to the kicks. It is worthly noticing that in the limit 
of e — > 0, Eq.© coincides with the result given by the 
Zurek's model. 

Even though this is a very simple quantum open sys- 
tem model, it is of great interest due to the fact that this 
scheme enables simulation of the quantum decoherence 
that usually appears for larger environments. As we have 
mentioned, one qubit as environment is not enough to 
produce decoherence on the system qubit in the Zurek's 
model. However, in the present case, the phase damp- 
ing is induced by a sequence of kicks that affect only the 
environment qubit, generated by a magnetic field that 
rotates the environment spin by an angle e. We believe 
that this practical implementation could be suitable for 
measuring of the complete GP in the case of a nonunitary 
evolution. 

This work was supported by UBA, CONICET, and 
ANPCyT, Argentina. 
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